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In this paper, we present a functional fluid limit theorem and 
a functional central limit theorem for a queue with an infinity of 
servers M/GI/co. The system is represented by a point-measure val- 
ued process keeping track of the remaining processing times of the 
customers in service. The convergence in law of a sequence of such 
processes after rescaling is proved by compactness-uniqueness meth- 
ods, and the deterministic fluid limit is the solution of an integrated 
equation in the space S' of tempered distributions. We then establish 
the corresponding central limit theorem, that is, the approximation 
of the normalized error process by a iS'-valued diffusion. We apply 
these results to provide fluid limits and diffusion approximations for 
some performance processes. 

1. Introduction. The queues with an infinite reservoir of servers are clas- 
sical models in queueing theory. In such cases, all the customers are imme- 
diately taken care of upon arrival, and spend in the system a sojourn time 
equal to their service time. Beyond its interest to represent telecommunica- 
tion networks or computers architectures in which the number of resources is 
extremely large, this model (commonly referred to as pure delay queue) has 
been often used for comparison to other ones whose dynamics are formally 
much more complicated, but close in some sense. Then the performances of 
the pure delay system may give good estimators, or bounds, of that of the 
other system. 

The studies proposed in the literature mainly focused on classical descrip- 
tors, such as the length of the queue: among others, its stationary regime 
under Markovian assumptions ([11]), the transient behavior and law of hit- 
ting times of given levels ([14]), the fluid limit and diffusion approximations 
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of normalized sequences ([3, 15]) are now classical results. Let us also men- 
tion the recent study on the existence/uniqueness of a stationary regime for 
the process counting the largest remaining processing time of a customer in 
the system ([21]). 

But when one aim is to have more accurate information on the state of 
the system, such as the total amount of work at current time {workload) or 
the number of customers having remaining processing time in a given range, 
no such simple state descriptor can be used. In order to address such ques- 
tions, one has to know at current time the exhaustive collection of residual 
processing times of all the customers present in the system. Consequently, 
we represent the queue by a point measure-valued process {fJ,t)t>Oi putting 
Dirac measures at all residual processing times. The price to pay to have 
such a global information is therefore to work on a very big state space, in 
fact, of infinite dimension. In the past fifteen years an increasing interest has 
been dedicated to the study of such measure-valued Markov processes (for 
definition and main properties, see the reference survey of Dawson [5] on this 
subject). Such a framework is particularly adequate to describe particles or 
branching systems (see [5, 22, 24]), or queueing systems whose dynamics are 
too complex to be carried on with simple finite-dimensional processes: the 
processor sharing queue (see [12, 26]), queues with deadlines (see [6] for a 
queue under the earliest deadline first service discipline without reneging, 
[7, 8] for the same system with reneging and [13] for a processor sharing 
queue with reneging), or the Shortest Remaining Processing time queue 
([1, 20]). 

In this paper we aim to identify the "mean behavior" of the measure 
valued process {iJ,t)t>o describing the pure delay system, introduced above. In 
that purpose, we use the recent tools of normalization of processes to identify 
the fluid limit of the process, or formal law of large numbers. This fluid limit 
is the continuous and deterministic limit in law of a normalized sequence of 
these processes. We characterize as well the accuracy of this approximation 
by providing the corresponding functional central limit theorem, that is, 
the convergence in law of the normalized process of difference between the 
normalized process and its fluid limit to a diffusion. 

Formally, it is rather straightforward in our case to identify the infinitesi- 
mal generator of the Feller process {iJ,t)t>o- The natural but unusual term is 
that due to the continuous decreasing of the residual processing times at unit 
rate as time goes on. This term involves a "spatial derivative" of the measure 
Ht, a notion which can only be rigorously defined within the framework of 
distributions. Because of this term, the fluid limit equation [see (4)] is the 
integrated version of a partial differential equation rather than an ordinary 
differential equation, as it is the rule in the previously studied queueing sys- 
tems. Thus, the classical Gronwall's lemma is of no use here. Fortunately, 
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we can circumvent this difficulty by solving the involved integrated equa- 
tion, known as transport equation, which is simple enough to have a closed 
form solution — see Theorem 1. Then, we can proceed using more classical 
techniques to show the convergence in law of the normalized sequence to the 
fluid limit (the authors have been informed during the review process of this 
paper that a similar result has been announced independently in [13]). As 
a second step, we show the weak convergence of the normalized sequence of 
deviation to the limit to a diffusion process. 

This paper is organized as follows. After some preliminaries in Section 2, 
we define properly the profile process {fJ-t)t>o in Section 3, show, in partic- 
ular, that (//f)j>o is Feller-Dynkin, and give the corresponding martingale 
property. In Section 4, we give the fluid limit of (/it)t>o- We deduce from this 
result fluid approximations of some performance processes in Section 5. We 
prove the functional central limit theorem for (//t)j>o in Section 6, and give 
the diffusion approximations of the performance processes in Section 6.3. 

2. Preliminaries. We denote by Dfy (resp. Ci,, Ck) the set of real-valued 
functions defined on R which are bounded, right-continuous with left-limit 
(rcll for short) (resp. bounded continuous, continuous with compact sup- 
port). The space Pfe is equipped with the Skorokhod topology and Cb with 
the topology of the uniform convergence. The space of bounded (resp. with 
compact support) differentiable functions from M to itself is denoted by 
(resp. C\) and for (p gC^, 

||0||oo:=sup(|0(x)| + |0'(x)|). 

We denote by 1 the real function constantly equal to 1, /, the identity 
function I{x) = x, x S R, and for all Borel set B G 53 (R), 1b the indicator 
function of B. For all f £T>i, and all x G R, we denote by Txf, the function 
rj{-):=fi--x). 

The Schwartz space, denoted by 5, is the space of infinitely differentiable 
functions, equipped with the topology defined by the semi- norms 

\(t>\f3,-f :=sup 

Its topological dual, the space of tempered distributions, is denoted by 
S' , and the duality product is classically denoted {fi,4>)- The distributional 
derivative of some fj, G S' is fi' £ S' such that {fi', cp) = —{fi, (f)') for all (p £ S. 
For all fi £ S' and x G R, we denote by Txfi the tempered distribution satis- 
fying {Txfi, <p) = {fJ-, Tx(t)) for ah (p^S. 

The set of finite nonnegative measures on R is denoted by A^^ (which is 
part of S') and Mp is the set of finite counting measures on R. The space 
A^t is equipped with the weak topology (T(A^t,Cfe), for which A^t is Polish 
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(we write = //(i^ for G and / G P;,). We say that a sequence 

{/i"}„gN* of J converges weakly to /i, and denote =4- /i, if for all / G C;,, 
(/x",/) tends to (//,/). 

We also denote for all x G M and all G A^j , r^;!/ the measure satisfying for 
all Borel set B, Txv{B) := v{B — x). Remark that these last two definitions 
coincide with that in S': for all (j)^ S' and v G ■M.'j , {v, (l))si = {u, (t>) , and 

for all X G M, {txI^, 4>) = {v-, T^cj)). 

For a random variable (r.v. for short) X defined on a fixed probability 
space (0,^), we say that converges in distribution to X, and 

denote ^ X, if the sequence of the distributions of the X^'s tends weakly 
to that of X. 

Let C(7Wj,IR) be the set of continuous functions from A^j to M. Let- 
ting < T < CO, for ii^ a Polish space, we denote C([0,r],i?), respectively 
P([0,r],£'), the Polish space (for its usual strong topology) of continu- 
ous, respectively rcll, functions from [0,r] to E. The mutual variation of 
two local martingales (Mt)t>o and {Nt)t>o in ^([0,T],M) is denoted by 
(< M,N >t)t>o with (< M >t)t>o = (< M,M >t)t>o for the quadratic vari- 
ation (or the increasing process) of {Mt)t>o- 

Since S' is a nuclear Frechet space (see [27]), note that {nt)t>o £ C([0, T],S') 
(resp. V{[0,T],S')) if and only if ((^t, (/'))t>o G C([6;r],M) (resp. 
V{[0,T],m)) for ah 0g5. 

Let T > 0, and let (Xt)t>o^ be a process of P([0,r],5'). We denote for 
all t G [0, T], /o ds, the element of S' such that, for all ^ G 5, 

l^j\,ds,(^^= j\xs,(k)ds. 

Let us, moreover, define the following 5'-valued processes: 

(1) {G{X)t)t>^ = {TtXt)t>o, 

(2) {H{X)t)t>o=( f'xsds) , 

(3) {N{X)t)t>o=( f\t-s{Xsyds] , 
which means, for all G 5 and all t G [0,T], 

{N{X)t,4>) = - f\Xs,Tt-s4>')ds. 

Jo 

Let {J='t)t>o be a filtration on (M,<B(M)). Let us recall (cf. [17, 27]) that a 
process M G P([0,T],5') is a 5' valued Tt-semi-martingale {resp. local mar- 
tingale, martingale) if, for all <j) £ S, {{Mt,(j)))t>o is a real ^^-semi-martingale 



^when no ambiguity is possible, we often write for instance X for the process {Xt)t 
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(resp. local martingale, martingale). According to [18], page 13, the tensor- 
quadratic process << M >> of the S' valued martingale M is given for all 
t > and all by 

« M V) :=< (Al, (t^),{M.,^) >t. 

For all 5'-valued semi-martingale M, and all real semi-martingale X, the 
5'-valued stochastic integral of M with respect to X is denoted by /g Mg dXg , 
and is such that, for all t > and all (j) £ S, 

''\MsA)dXs. 









)-L 







Let T > 0. We say that X G ^([O, T], 5') satisfies the integrated transport 
equation associated to {K,g) (see [8]) if, for some K (z S' and {gt)t>o £ 

(4) Xt = K+ f\xs)' ds + gt for all t G [0, T] . 

Jo 

Let us then recall the following result. 

Theorem 1 ([8], Theorem 1). The only solution in I){[0,T],S') of the 
integrated transport equation associated to {K,g) is given for all t G [0,T] by 

(5) Xt = TtK + gt + N{g)t, 
where the mapping N is defined in (3). 



Note that Propositions 3 and 4 show, in particular, that X defined by (5) 
is an element of T>{[0,T],S'). 

3. The profile process of the M/GI/oo queue. Hereafter we consider a 
pure delay system M/GI/co: on a probability space {^l,T,P), consider a 
Poisson process 

Nt ■■= J2 ^{T^<t} 

iGN* 

of positive intensity A, representing the arrivals of the customers in a queue- 
ing system with an infinite reservoir of servers. Hence, any of them is im- 
mediately attended upon arrival. The time spent in the system by the ith 
arriving customer (denoted Cj) equals the service duration fjj he requests. 
We assume that the sequence of marks {uijigN* is i.i.d. with the random vari- 
able a having the distribution a G ■M.'f ■ We denote at all time t>0,Xt the 
number of customers currently in the system, and St the number of already 
served customers at t, related to Nt and Xt by the relation Nt = Xt + St- 
The workload process {Wt)t>o equals at time t the total amount of service 
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requested by the customers in the system at t, in time units. The profile pro- 
cess of the queue is the point-measure valued process {nt)t>o whose units 
of mass represent the remaining processing times (i.e., time to the service 
completion) of all the customers who already entered the system at current 
time. In other words, for all t > 0, 

Nt 

^* = X!%+'^^-^)• 
^=l 

In this expression a nonpositive remaining processing time (Tj + — t < 0) 
stands for a customer who already left the system, precisely t — (Tj + cjj) 
time units before t. As easily seen, for all t > and any x <y, the number of 
customers having at i a remaining processing time in {x,y) (the boundaries 
may be included or infinite) is given by {fit,'^{x,y)) ■ Thus, the congestion, 
service and workload processes X, S and W are easily recovered from fi by 
the relations Xt = {fj-t, Ir^), St = {fit, 1r-) and Wt = {fj-t, I1r*_^) , t > 0. 

The dynamics of {nt)t>o can be described as follows. For any t > 0, fj-t+dt 
is the translated of fit toward left by dt, that is, T^tfit (the remaining service 
times of all the customers in the system at t, if any, decrease at unit rate), 
and a Dirac measure at a is added if a new customer enters the system 
between t and t-\-dt having service time a. In other words, we have for any 

(6) {f^t+dtA) = {Tdm,(l)) + (l){a)N{dt,da), 

where N{-,-) denotes the Poisson measure associated to the marked arrival 
point process. Hence, the above equation allows one to construct {{nf,4>))t>Q 
for any (f) G Ck by induction on the arrival times and for any initial state 
fj, £ A4~j (as in [26], page 192). Since Ck is a separating class of TWj , the 
initial value /xq and the evolution equation (6) fully define {iJ,t)t>o- Let h> 
and F be a bounded continuous function from Ai^ into M. Denote for all 

Ai{h) = {There are exactly i arrivals in [0, h]}. 
We have for all fi£ M'j 

(7) =EE 



i>0 



■ V .7>0 / 



A,ih) I ^0 - M 



(1 - \h)F{Thii) + \hj F{Thifi + 4)) da{x) + e(/u, h), 
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where 

e{fi,h) = {P[Ao{h)] - (1 - Xh)}F{Thf,) 

+ {P[Aiih)] - Xh} J F{Th{fi + 6,))da{x) 



(8) 



i>2 



Hence, e(/i,/i) is a o(/i), as easily seen using dominated convergence for the 
third term on the right-hand side of (8), and from the fact that F is bounded 
for the other three terms. Therefore, according to [5], page 18, {iit)t>o is a 
weak homogeneous A^j -valued Markov process having transition ThF for 
all h > and all bounded continuous F. Note, moreover, that {fJ.t)t>o £ 
V{[0,oo),Mp since ((/xt, 0))t>o G P([0, cx)),M) for all cl)£Cb. 

Proposition 1. The process {fJ-t)t>o a Feller-Dynkin process o/'D([0, 
oo),M+). 

Proof. According to Lemma 3.5.1 and Corollary 3.5.2 of [5], the Markov 
process {fit)t>o enjoys the Feller-Dynkin property if: 

(i) For all / G Cj^, h> 0, the mapping fi E[Fj(^/j) | hq = ^t] is contin- 



uous. 



(ii) For all h>0, E[Fi(/i/j) \ fio = fi] — >0, as /i(M) +oo. 
(in) For all /j, G M~f and / G C}^, E[Fj(^/,) | /J-o = — > Ff{fi) as h^O. 

Let us denote for all / G Cj^, Ff the mapping from into M defined by 
view of (7) and (8), we have for all /i G Ai~j and / G Cj^ 

B[Ff{fXh)\fio = fA 



(9) 



where 



ef{h) = {P[Ao{h)] - (1 - Xh)} + {P[Ai{h)] - Xh} J e 
+ E[{e-/(--('^-^^))-e-/(--^)}U(,)] 



-f{x-h) 



da{x) 



i>2 L Kj=l 



n 
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which is a o{h) from the same arguments as for (8). Note, moreover, that 
€f{h) does not depend on Hence, the continuity in (i) is granted by the 
fact that, for all h> 0, the mapping ;U i— > t/j^ is continuous from into 
itself. For (ii), remark that the total mass on R of ThfJ- equals that of /i. 
Finally, the convergence in (iii) holds since {fi,Thf) — >h^o (/^j/) by domi- 
nated convergence. The proposition is proved. □ 

As easily seen from (9), and since the set of linear combinations of Fj, 
f G C]^, is dense in C{M-~j ,M.) (see [26], Proposition 7.10), the infinitesimal 
generator A of {fJ.t)t>o is given for all fi G by 

h--*0 h 
= lim ^(■'^hfi) -F{^i) _ r ^ ^ 

h~*0 h J 

for all F S C(A^j , M) such that the latter limit exists (we say that F belongs 
to the domain of A). 

Proposition 2. For all </> G C^, the process defined for all t>0 by 
Mt{4>) = {ixt,4>) - il-ioA) + I {UsA') ds - Xt{a, 4>) 



is an rcll square integrable J^t-n^'iirtingale. For all ^S, the mutual vari- 
ation of (Mt{(p))t>o with {Mt{ip))t>Q is given for all t>0 by 

(10) <MX(l)),M.i^P)>t = Xt{a,cj)iP). 

Proof. For all (/> G C^, the mapping II^ from M'j into M defined by 
n0(/x) = clearly belongs to the domain of A. As a consequence of 

Dynkin's lemma ([9, 10]), the process defined for all t > by 

Mt((/.)=n^(/it)-n^(/io)- fAU^{fis)ds 

(11) 

= {lJ't,(f>) - {fJ-o, <P) + {lJ's,(t>') ds - Xt{a, (j)) 



is an rcll ^t-local martingale. Let now ijj ^C^. The mapping II,^^^ from A^j 
into M defined by n^^^(/i) = also belongs to the domain of A, 

implying that 

(12) Mt((/.,V') = n<^(Ait)n^(Ait)-n0(/xo)n^(;uo)- / Au^,^{fis)ds 

Jo 
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is as well an rcll JT^-local martingale. But as easily checked, for all /x G > 

(13) AU^,^{fi) = U^{n)An^{fi) + U^{fi)AU^{fi) + X{a, (1)^}). 
Ito's formula yields together with (11) that, for all t > 0, 

n<^(//t)n^(/it) 

= n0(/io)n^(/uo) + / n0(/is)dMs('0) + / ii^{ixs)Aii^,{ixs)ds 

Jo Jo 

+ f n^(/i,)dM, ((/.)+ / n^{fis)Au^{fis) ds 

Jo Jo 

+ <n^(Ai.),n^(/i.) >t, 

which, together with (12) and (13), implies 

/ u^{ns)dMsW+ f u^{^ls)dMs{(|)) + <u^i^i.),'^i>{^^.)>t 

Jo Jo 

By identifying the finite variation processes, we obtain that P-a.s. for all 

t >0, 

(14) <n0(;u.),n^(/i.) >t = At(a, (/)?/'), 

but in view of (11), this last quantity equals < M.(i;^)), ^^(-0) >t. In partic- 
ular, for all t > 0, 

E[<M.((/)) >i] =At(a,(/)2) <cx), 
hence, (Mj(</)))t>o is a square integrable martingale. □ 

4. Fluid limit. Consider a sequence of M/GI/oo systems such that the 
nth system has an initial profile /^q, is fed by a Poisson process {Np)t>o of 
arrival times {T"}jgpj» and of intensity A", in which the customers request 
i.i.d. service durations {c"}jgN* which have the nonatomic distribution a" of 
a r.v. 0"". We assume, furthermore, that o"'" is integrable, that is, (a", I) < oo. 
The process {fif)t>o, defined for all t by 

i=l 

is the profile process of this nth system. Denote by (^")t>o the associated 
filtration. Let us also define as previously the performance processes of the 
nth system: X", 5" and W^. We normalize the process in time, space 
and weight by defining for all t > 



10 



L. DECREUSEFOND AND P. MOYAL 



Thus, for all Borel set B and all t, 

n 

where nB := {nx, x G B}. In words, p,^ isl/n times the point measure having 
atoms at levels n times smaller than that of //J^j. We also define (^p)j>o := 
{J-^f.)t>o, the associated filtration, and normalize the arrival process as well 
as the performance processes of the nth system the corresponding way, that 
is, for all t > 0, 

j\jn _ vn 

n n + 

_ on _ Ti/n 

5r := ^ = (Hi, 1m-), ■.= -f = ifl^J^Mx)- 

n + 

We denote for all i G N*, := /n. Thus, the {a^jigN* are i.i.d. with the 
nonatomic distribution a" of the r.v. a'^/n. The distribution a" is such that 
{a^,I) < oo, and satisfies for all Borel set B 

a'^{B) = a''{nB). 

We assume hereafter that the following hypothesis holds. 

Hypothesis 1. There exists A > such that 

(15) A" A, 

n— +00 

For all e > 0, there exists > such that, for all ?i G N*, 

(16) P[(^f^,l)>nM,]<e. 
For some measure JIq oi M.'^ such that 

(/io,/) < oo, 

for ah / G Pfe, 

(17) i^jinj) —,Jj,lJ) in probability. 

There exists a nonatomic probability distribution a* such that 

a" ^ a*, 
(a",/) — > (a*, I) <oo. 

Let 4>^Cl and '0"(") = (j^i- /n)/n. According to Proposition 2, the process 
defined for all t by 

(18) Mr(</'):=M:,(V'") = (Ar,</')-(/i^,<A)+ r(/^-'^Vs-A"t(a",<A) 
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is a square integrable ^"-martingale of ^([O, oo),M), such that 
(19) <M"((A) = 

n 

In other words, the process M" defined for all t by 





is a 5'- valued ^"-martingale of tensor-quadratic process given for any </> and 
in 5 by 

«M"»t ((/),V) = — t(«",</'V')- 
n 



Theorem 2. Assume that Hypothesis 1 holds. Then 
Jf'^ji* mP([0,oo),7Wj:), 

where p.* is the deterministic element of C{[0,oo),M^) defined for aUt>0 
and all <j) GDi, by 



(20) {fllcj)) = {fllTt<l))+X f\a*, 

Jo 



t 

, .Ts(l))ds. 





Proof. First we prove that {/i"}neN* is tight in 2?([0,oo), A^jr). To 
this end, it suffices to show conditions C.l and C.2 of Theorem A. 2 in the 
Appendix. To show C.l, fix <^ G and T > 0. Remarking that, for all s > 0, 

equation (18) yields P-a.s. for all u <v <T, 

< <p')\ds + A"|(a", </.)||7; -u\ + \M;{v) - M^{u)\ 

(21) ^" , _ 

<b-n|{||0'|U(iV? + (Ao,l)) + ll'/'l|ooA"} 

+ \M;{v)-M^{u)\. 
Let ^ > 0. From Doob's inequality. 



P 



.t<T 



n—roo 



in view of (19). Hence, it follows from the standard convergence criterion 
on P([0,T],]R) that {(M"(0))f>o}neN* converges in distribution to the null 
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process. This sequence is, in particular, tight, and it is routine in view of (21) 
and Hypothesis 1 to check the standard tightness criterion in P([0,T],]R): 
for all e > and i] > 0, there exists 5 > and N G'N such that, for all n > A^, 



(22) 



sup \{fi:„c^)-{fi':„cP)\>r^ 
■u,v<T,\v—u\<S 



< £. 



Letting (3e ■ = M^\\(j)\\oo, assumption (16) implies that, for all n G N*, 
(23) Pmo^^)\ > Pe] < P[||</'l|oo(/iS', 1) > M,U\U < e. 



Hence, (22) and (23) show that {{{flt,(t>))t>o}n€N' is tight in V{[0,T],R) for 
all T > (see [25], Theorem D.9). This sequence is thus tight in V^O, oo),M). 

We now prove condition C.2 (compact containment). Fix T > 0. Let us 
first apply [12], Lemma A. 2.: under Hypothesis 1, we have 



nNV 



(24) \[-T.<^i-^) 



iXt{a*,<P))t>o inP([0,r],M) 



i=l 



t>0 ) nSN* 



for any measurable cf), such that (f) is continuous on the supports of a* and a", 
n E N*, and such that {\(j)\,a*) < oo and (|(/)|, a") < oo, n G N*. In particular, 
this yields for any < Z < T, and any such (p, 



(25) 



1 



N" 



sup - </'(^»") >2A/(a*,( 

te[o,T-/] n ^ 



0. 



Taking I = T and = 1 (resp. (f) = I) in the above expression yields 

Pim > 2XT] — > 0, 

n— >cxD 



1 



Y^a?>2\T{a*,I) 



n " , 



0. 



Letting Mt = ma^{2 XT + (fi^,!), 2 XT {a* , I) + (/x^,/)} + 1, we have 



sup max{{fil\lM+),{fiUlK+)}> Mt 
■te[o,T] 



(26) < P[(Ao", 1) > {flo, 1) + 1] + mr > 2Ar] 



-Y^-a->2XT{aM) 

i=l 



0. 



n— >+oo 



For all < 7? < 1 , denote 

/Ct,^ := {C G A^;;max{(C, 1k+), (C,I1r+)} < Mt, (C, 1(-oo,-T]) = 0}. 
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Since (C,/1m+) < Mt implies that for all y > 0, Cibli oo)) < Mx/y, we have 
lim sup C([y,oo)) = 0, lim sup (((—oo,y]) = 0, 

which implies in turns that KLtti C M.~j is relatively compact (see [16]). Now, 
since up to time T no already served customer can have a remaining pro- 
cessing time less than — T, we have 

SUp(p",l(_oo -T]) =0. 
t<T 

This, together with (26), implies that 

liminf P[/ir G /Ct«, for all t £ [0,T]] > 1 - r?. 

Kt,j; being the closure of /Ct,?;, we found a compact subset KT,ri such 
that 

liminf P[/xr G Kt„, for all t G [0,rl] >l-r], 

which completes the proof of tightness. 

Let now (xt)t>o be a subsequential limit of {/x"}nGN*- We have for all 
t > 



0+ f iXs)' ds + Xta* 
Jo 



Xt = IJ-o 

which is nothing but the integrated transport equation associated to {flQ,g), 
where gt ■= Xta* , t > 0. From Theorem 1, we have for all t > and (j) £ S 

{Xt,(p) = {nflo, 0) + {9t,(l>) + {Nigt),(p) 

integrating by parts. The subsequential limit is therefore unique in the space 
^([O, oo), ), and equal to /i*, since 5 is a separating class of -M^ . □ 

5. Fluid limits of some performance processes. Let us provide some ap- 
plications of Theorem 2 to the asymptotic estimation of some performance 
processes describing the queueing system. Assume that Hypothesis 1 hold 
and, in addition, that the limiting initial profile fi^ in Hypothesis 1 is such 
that 

(27) has no atom. 

Assumption (27) implies in view of (20) that fi^ has no atom for any t > 0. 
Thus, Theorem 2 and the Continuous Mapping Theorem yield that, for 
any x <y, the sequence {((/^f, l(x,y)))t>o}r!,eN* tends in distribution to the 
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real deterministic function ((/ij', Ij-^. y)))t>o (the boundaries may as well be 
included or infinite). Thus, a fluid limit approximation of the process count- 
ing the customers having residual service times in the range {x,y) is given 
for alH > by 

{fJ'h'^{x,y)) = ^''oii^ + t,y + t)) + XJ^ l^J ^ da*{x)jds. 

In particular, a fluid approximation of the normalized congestion process 
X" is given by the process X* defined for all t by 

X; = ifll = fioiit,^)) + Xj^ da*{x)^ ds, 

whereas the normalized service process S"" can be approximated by S* , 
where 

5; = {fil, 1k_) = /i5((-oo,t]) + ^^*(^' da*{x)^ ds. 
Remark now that for all t <T, all x > and all n G N*, 

1 nt 

wr < (Ar,^i(o,x]) + (Ao,^i(i+x,oo)) + -E^ri^r>- 

1=1 

In view of (24), the third process on the right tends to (At(Q;*, Jl(^. oo)))t>0- 
Hence, it is bounded over compact sets and uniformly in n since (a*, 00)) 
is finite. The same argument applies to the second process on the right 
since (/Iq , is finite as well, and the first one tends in distribution 
to ((/Uj , /l(o,a;]))t>o- Heuce, {W^^lnGN* is tight, and any subsequential limit 
thus reads ((/ij , /1k;^))j>o. In other words, the fiuid limit of W"^ is given by 
W* , where 

n+00 ft 



{x — t) dfiQ{x) + X J {x — s) da* (x)^ ds . 



6. Functional central limit theorem. 



6.1. Preliminary results. In this section we prove two technical results 
(Propositions 3 and 4), which will be useful in the sequel. We refer again the 
reader to the definitions and notation introduced in Section 2. Throughout 
this whole section, fix T > 0. 

Lemma 1. For all (p gS, the collection {Tr4>,r £ [0, T]} is a pre-compact 
set of S . 
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Proof. For all semi-norm associated to the integers /?, 7, for all r S 

\rr<i>\p,, = sup|x^</.W(x - r)| < ^ Cy-^\<l>\,., < ^ C^T'^"^ |</>|,,^. 

i=i i=i 

The set {Tr4>,r G [0, T]} is thus bounded in the nuclear space 5: it is a 
pre-compact set in view of Proposition 4.4.7, page 81 of [23]. □ 

Lemma 2. For all </> G 5, the mapping r \ — > Tj.^ is continuous from M 
into S. 

Proof. Let r' < r. Let /3 and 7 G N. For some s G ]r',r[, 

\Tr'(l) — Tr4>\i3^-y = SUp|x^((/)^'^^ (x — r') — 4>^'^\x — r))\ 

X 

<\r-r'\j2cy\4>\,,,+, 
i=i 

P 

i=i 

Thus, for all e > 0, fixing := ^(Mr^^^^^^)"^, we have \Tr'4> — Tr(t>\f3^'y < £ 
whenever \ r — r' \ < r]. □ 

We therefore have the following results. 

Proposition 3. The mapping G defined by (1) is continuous from 
C{[0,T],S') into itself 

Proof. We first prove that G takes its values in C([0, T],S'). Let </> G 5, 
X eC{[0,T],S'), and let G^ be the mapping from into M defined by 
G^{t,r) = {Xt,Tr(t)). Let t,r<T. For all t',r', 

(28) \G^{t',r') - G^{t,r)\ < \{Xt, - Xt,Tr^)\ + ^c/. - r,(/>)|. 

On the one hand, {{Xs,(p))s>o G ^([0, T], M), and is thus bounded on [0,r]: 

sup|(Xs,(/>)| = sup \{Xs,(p)\ < oo. 
s<T se[o,T]nQ 

Therefore, X := {Xs,s G [0, T] n Q] is a weakly bounded subset of the set of 
continuous mappings from S into M. The space S is Frechet, and hence is a 
tonnel in view of corollary 0, page IIL25 of [4]. From the Banach-Steinhaus 
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theorem (Theorem 1, page III. 25 of [4]), X is equicontinuous. Hence, for all 
e>0, 

Ve:= n ^s'(]-'7;,'J) = \'P, sup \{Xs,cP)\<^\ 



sG[0,T]nQ 



is a neighborhood of 0, the zero of S. The mapping r t,.(/) being continuous 
in view of Lemma 2, for some rj^^^ > 

\r' -r\< rjr^^ =^ {Tr'(j) - Tr(j)) G K =^ \ {Xt,Tr'(j) " Trt/')! < |. 

On the other hand, {{Xs,Tr(l)))s>o G C([0,r],M), hence, for some Sr^t,<f> and 
all t', 

\t' -t\ <5r,t,<P^\{Xt' - Xt,Tr4>)\ < |. 

For all t, r and all e > 0, the two previous relations in (28) imply that, for 
some ii^r,4, and 5r,t,(t>, 

\t' -t\< 6r,t,^, \r' -r\< rjr,,f> =^ \G^{t',r') - G^(t,r)\ < e. 

The mapping is thus continuous, and hence uniformly continuous on the 
compact set ([0,T])^: for some 5^,ri^>Q, 

sup sup \G^{t',r')-G^{t,r)\<e. 

r,r',\r-r'\<ri^ t,t'\t-t'\<S^ 

Finally, letting C</, := r?^ A 5^, 

sup \{G{X)t>,<P)-{G{X)t,(t>)\ 

t,t' ,\t-t'\<i^ 

< sup sup \G^{t\r')-G^{t,r)\<e. 

r,r',\r-r'\<r]^ t,t' ,\t-t'\<5^ 

Thus, G{X) £ C{[0,T],S') since the map {G{X)t.,4>) is continuous for all 

Let us now show the continuity of G. Let be a sequence of 

P([0,r],5') tending to X. In particular, for ah (/> G 5, {((X^", (/)))t>o}nGN* 
tends in ^([O, T], M) to {{Xt,(p))t>o, and thus sup^^j^. supf<^j:\{Xp,(p)\ <oo. 
Hence, in view of the Banach-Steinhaus theorem, 

7t::={Xf,nGN,tGQn[0,r]} 

is equicontinuous. Therefore, for all (j) £ S, 

snp\{X^-Xt,<P)\ -^0, 
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and consequently, for all k, all collection = 1, . . . , A;} of 5 and all i <k, 
sup max \{XV- - Xt, (j)^] — > 0. 

i=l,...,fc n^oo 

The latter means that for all semi-norm of the weak topology, 

snppUX^-Xt) 0. 

t<T ' n^oo 

Thus, in view of Pietsch's theorem ([23], Proposition 0.6.7, page 9), for all 
pre-compact set K of 5, 

supsup|(X,"-Xt,(/))| 0. 

This result can be applied for a fixed (/> G 5 to the collection {T,.(/),r G [0,T]} 
(which is a pre-compact set of S from Lemma 1) to obtain 

supsup|(X4"-Xt,r,(/))| 0. 

t<Tr<T 

It follows that G is continuous, since 

sup I - 0) I = sup I - Xi, rt<A) I 

t<T t<T 

< sup sup I {X^ - Xt,Tr(l))\ ^ 0. ^ 
t<Tr<T n^co U 

Proposition 4. T/ie mapping N defined by (3) is continuous from 
V{[0,T],S') into C{[0,T],S'). 

Proof. That N takes values in C{[0,T],S') can be shown similarly to 
Proposition 3. For the continuity of A'^, remark that the mapping H defined 
by (2) is continuous from I?([0,T],5') into C{[0,T],S') since, as well known, 
the mapping {Xt)t>o ' — > {Jq Xs ds)t>o is continuous from 'D([0,r],M) into 
C([0,T],M) for the Skorokhod topology. The proof then proceeds as that of 
Proposition 3. □ 

6.2. Central limit theorem. Fix {A"}„gN* ) A, and {cr^lngj^* satisfying 
Hypothesis 1. {p-"'}neN* is the corresponding sequence of normalized pro- 
file processes. From Theorem 2, tends in distribution to jl* in 
T){[0,oo),A4~j), defined by (20). Throughout this section, we make the fol- 
lowing additional assumptions. 

Hypothesis 2. 
(29) V^(A" - A) — ^ 0. 

n — ^oo 
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There exists yo £ S' , such that for all (j) £ S, 

(30) Vn{jlo -K^4>) — > (3^0, in probability, 

(31) v^|(a",/) -(a*,/)| ^ 0. 



Let n G N* . Define the following processes for all t by 

^r(<^,V') = A"t(a",(/>V^), <^,Ve5. 

The process A^" is an 5'-valued ^"-martingale, hence, for all (f) £ S, the 
process {{J^f , (p))t>o is a real square integrable martingale of increasing 
process (^J^((/., (/.))t>o [by (19)]. 

Lemma 3. Under Hypotheses 1 and 2, for all T > 0, the sequence {A^'"}„gN* 
converges in distribution to the martingale Ai o/D([0, T], 5'), whose tensor- 
quadratic process satisfies for all t £ [0, T] and all £S 



(32) 



«M >>f((/),V') = Xt{a*,(j)'il'). 



Proof. We use the convergence criterion, Theorem A. 3 in the Appendix. 
Define for all t > and (j),ip £S 

7t('A,V') := At(a*,(/)V'). 

For ah G 5, 7o(0, 0) = and for ah n, (((7W^, 0))2 - (</., cp)\>Q is a ^f- 
local martingale, thus conditions (37) and (38) are verified. On the other 
hand, for all T > and ip £S, we have P-a.s. for all t < T 



A'}{(t>, V) - iM, V') = (A" - ma^'^H) + At(a" - a*, </>V'>, 



hence, 



supMr('A,V')-7t('/',V')}' 



< 2r2{(A" - A)^ 



+ A^ 



n\2 
lloo 



(a" - a*,/)} 



0, 



which proves (41) taking ip := (p. Finally, for all T > 0, condition (40) is met 
since {A2{(pA))t>o C{[0,T],R) and 



E 



sup((A^r,</')-(A^L,'A))' 

t<T 



< 



E 



sup((/ir,<A)-(Ar-, </>))' 



n n— >oo 
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the jumps of {{[If , (j)))t>o being a.s. of size less than Mh^- have (39) as 
weh. Therefore, we can apply Theorem A. 3 to state that 

(33) {{{M^,cl)))t>o}neN'^{Pt)t>o inP([0,cx)),M), 

where P is a continuous martingale of increasing process (7j((^, (/>))t>o. In 
particular, for all T > and ah (?i G 5, {((A^?, </'))t>o}neN* is tight in V{[0, T],R), 
and hence (Theorem A.l in the Appendix), {M^}nm* is tight inP([0,T],5'). 
All subsequential limit M thus satisfies {Mt,<P) = Pt for all t G [0,T] and 
all G 5, which means with (33) that {A^"}„gN* tends in distribution to 
the martingale M of P([0,r],5'). Finally and according to the previous 
arguments, 

{(< (Air,(/.),(Air,V') >t)t>oW = {(^r(<^,V'))i>oW 

tends to (<< >>(((/>, V'))t>o as well as to (7j((/), i/;))t>o, which completes 
the proof. □ 

Consider now the Hilbert space L^((ia*). Since a* is a probability mea- 
sure, we have S C L^((ia*), and hence, (L^((ia*))' C S' . Moreover, L^((iQ*) 
admits a countable basis, denoted {^i}igN- For alH > and all if^ G L^((ia*), 
denote 

Ci('0)= j ip{x)hi{x)da*{x), 
the ith coordinate of tp in that basis. We have the following result. 

Theorem 3. Under Hypotheses 1 and 2, for all T > 0, the sequence 
{3^"}nGN* tends in distribution in P([0, r],5') to the process y* defined for 
all t>0 and all (pG S by 

(34) {y;,ct)) = {yo,Ttcl)) + VxY, f Ci{Tt-s^)dBi, 

where {-B*}i>o is a sequence of independent real standard Brownian motions. 

Proof. For all ?i G N and ah t G [0,T], P-a.s., 

y^ = y^ + [\y^y ds + y^{A"ta" - Xta*} + • 
Jo 

Hence, {yi^)t>o solves a transport equation, and in view of Theorem 1, the 
above equation amounts P-a.s. for all t < T to 

y^ = G{y^)t + y^{A"ta" - Xta*} + TW^ 

+ V^A"iV(- a") - V^\N{- a*) + N{M'')t 

(35) = G{y^)t + V^(A" - A)ta" - Xt^/E{a* - a") + M'l 

+ Vn(A" - X)N{- a'')t - X{N{^/^ ■ a*)t - ■ a")t} 

+ A^(A^")t, 
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where (3^o )t>o denotes the 5'-valued process constantly equal to 3^q and • a" 
(resp. -a*, -^n-a", ^/n • a*) denotes the ^'-valued process {ta"')t>o [resp. 
{ta*)t>o, {y/ntci^)t>Q, {y/nta*)t>Q\- Fh'st, from (30), for all S 5, ^ > 0, 

P[K3^o -3^o,<A)l>e]„-^0. 

Hence, in view of Corollary A.l in the Appendix, we have that (J'q )t>o ^ 
(3^o)t>o in ^([0,7"], 5'). Thus, from Proposition 3 and the Continuous Map- 
ping Theorem (see [2]), 

{G{y^)t)t>^ {Giyo)t)t>o in Vi[0,T],S'). 

On the other hand, we have 

sup V^|A"-A|t|(a",(/<)| < V^|A" -A|r||(/.|U ^0, 

0<t<T 



n— >oo 



sup V^|A" - A|KAf(-a")t,0')| < V^|A" - X\T'U'\\oo 0, 

sup AtV^|(a*,(/.) - (a",(/>)| < Arv^||(/.'||ooKa*,/) - ^0. 

0<t<T n^oo 

Hence, the sequences {(A" — A)-^/^ • a"}„gN* , {(A" — A)A^(v^- Q;")}„gN* and 
{A(-v/n-a* — \/^-a")}r!.GN* converge to the null process (0)(>o of C([0, T], 5'). 
The convergence of this last sequence implies in view of Proposition 4 that 

A{iV(V^.a*)-Ar(^.a")} ^ (0)t>o in C([0,T],5'). 

Now, from Lemma 3, =^ M in 'D{[0,T],S'), thus, in view of Proposi- 
tion 4 and the Continuous Mapping Theorem, N{M"-) =^ N{M) in I?([0, T] , 5') 
Consequently, {y^}neN* is tight and from (35), its limit in distribution y* 
satisfies P-a.s. for all t S [0,T]: 

(36) y; = G{yo)t + Mt + N{M)t = Ttyo + Mt+ fn^siMs)' ds. 

Jo 

In view of (32), the process M reads M = VXM{a*), where B(a*) is the 
cylindrical Brownian motion on (L^((iQ:*))' (see [18]). Thus, from (36), we 
have for all G 5 and for all t 

{yiA) = {yo,Tt(i}) + Vxl^{'B{(^it,ci>)- J\nais,Tt-s^')ds 

= {yo,rt(p) + V\y2{c,{(P)Bl - f\i{Tt^s(t>')Blds 

where {-B*}j>o is a sequence of independent real standard Brownian motions 
(see again [18]). Note that the latter series converges in L^(r2) since (j) £ 
Ij'^{da*). Remark now that, for all i > 0, 

d{ci{Tt^s(l))) , ,,^ 
ds = ^^(^*-'^)- 
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Using Ito's integration by parts formula, this yields to 

c^{^)Bl - f c,{Tt-s<t>')Bi ds = f Ci{Tt.s4>) dBl 
Jo Jo 

which completes the proof. □ 

6.3. Diffusion approximations of the performance processes. In this sec- 
tion we show that Theorem 3 can be used to provide diffusion approxima- 
tions for the congestion, service and workload processes. The stochastic in- 
tegrals on the right-hand side of (34) make sense for any (/> G L^(dQ:*), hence, 
the process (3^t*)t>o takes values in (L^((ia*))', provided that the limiting 
initial state 3^o belongs to (L^((ia*))'. Showing that in that case the con- 
vergence announced in Theorem 3 holds in P([0, T], (L^((iQ:*))') is an open 
problem at this point, that is beyond the scope of this paper. However, it is 
easily seen that under Hypothesis 1 the limiting service time distribution a* 
is such that the functions Ir^ and 1k_ belong to L^(da*). Hence, diffusion 
approximations for the congestion and service processes can be provided by 
the real processes obtained when fixing (j) = and (f) = 1r_ in (34), that 
is, for all t > 0, 

(3^;, 1m;) = yo{{t, +oo)) + VaJ^ ^* (^j^^ hi{x) da*{x)^ dBl 

and 

{ylU_)=yo{{-^,t]) + VxJ2 fit ' hi{x)da*{x)]dB\. 

Moreover, provided that f^°°x'^dQ*{x) < oo (i.e., the limiting service time 
has a finite second moment), a diffusion approximation for the workload 
process is given for all t > by 

p+oo 

(K,/lM;) = y^ {x-t)dyoix) 




Assume, for example, that a* is the exponential distribution e(l) [which is 
the limiting distribution obtained when taking cr" ~ e(l/n) for all n E N*]. 
In that case, it is well known that a possible {/ijjjgN is given by the sequence 
of Laguerre's polynomials, defined for alH > by 
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Then a diffusion approximation of the congestion process (and accordingly, 
of the service and workload processes) is given for all t > by 

APPENDIX: TIGHTNESS AND WEAK CONVERGENCE ON 
FUNCTIONAL SPACES 

Let us recall some results about tightness and weak convergence on metric 
spaces. 

Theorem A.l (cf. [19], Theorem 4.1). Let F be a nuclear Frechet space, 
F' be its topological dual, and {X^}n£fq* be a sequence of VdO^T], F') -valued 
r.v. Then {X"-}nm* tight if for all (p G F, the sequence {{{X't , (p))t>o}n&N* 
is tight inV{[0,T],R). 

This shows, in particular, since 5 is a nuclear Frechet space, that the tight- 
ness of {X^}n£N* in ^([0,r],5') is granted by that of {((X", (p))t>o}neN* for 
all S. Moreover, since 5 is naturally a separating class of its topological 
dual S' , we have the following: 

Corollary A.l. For all sequence {X'^}nm-' ofV{[0,T],S'),X''^X 
if for allcPGS, i{X^,cP))t>o ^ {{Xt,<P))t>o in V{[0,T],S') . 

The space has no such good topological properties as that of S' . 
Nevertheless, a tightness criterion is established for sequences of A^^ -valued 
processes. 

Theorem A. 2 (Jakubowski's criterion, cf. [5], Theorem 3.6.4). The se- 
quence {X^l^gN* of V{[0,T],A4^) is tight if the following two conditions 
hold: 

C.l. For all (j) ^ Cl , the sequence {{{X'^ ,({))) t>o}n€N* is tight in V{[0,oo),M.) , 
C.2. For all T > and <rj < 1, there exists a compact subset K^.r? of 
such that 

lim^inf P[Xf G Kt,,, Vt G [0, T]] > 1 - rj. 

Let us finally give the following result, which establishes a criterion for 
the convergence in distribution of a sequence of real processes to a diffusion 
process. 
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Theorem A. 3 ([10], Theorem 1.4). Let {i?"} and {A"} be two se- 
quences o/ P([0, oo), M), and an increasing function 7 G 2?([0, oo), M) such 
that 



(37) 



i?"(0) = and 7(0) = 0, 



(38) i?" and (i?")^ — A" are two Tt — local martingales, 

and for all T > 0, 



(39) 



E 
E 



sup{i?"(i) - i?"(i-)}^ 

t<T 



svip{A''{t)-A''{t-)} 

.t<T 



0, 



0, 



(40) 

(41) Vt>0,Ve>0 P[{^"(t) -7(t)} >e] — ^0 

n-~*oo 



Then, tends in distribution to a continuous martingale of increasing 
process {l{t))t>Q. 
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